Stress relaxation in a perfect nanocrystal by coherent ejection of lattice layers 
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We show that a small crystal trapped within a potential well and in contact with its own fluid, 
responds to large compressive stresses by a novel mechanism — the transfer of complete lattice 
layers across the solid- fluid interface. Further, when the solid is impacted by a momentum impulse 
set up in the fluid, a coherently ejected lattice layer carries away a definite quantity of energy and 
momentum, resulting in a sharp peak in the calculated phonon absorption spectrum. Apart from 
its relevance to studies of stability and failure of small sized solids, such coherent nanospallation 
may be used to make atomic wires or monolayer films. 

PACS numbers: 



Solids subject to large uniaxial deformations, re- 
lieve stress either by the generation and mobility of 
dislocations [l[ and/or by the nucleation and growth of 
cracks [H What is the nature of stress relaxation 
when conditions are arranged such that these conven- 
tional mechanisms are suppressed? Nano-indentation 
experiments Q show that if small system size prevents 
the generation of dislocations solids respond to ten- 
sile forces by shedding atoms from the surface layer. In 
this Letter, we study equilibrium and dynamical aspects 
of this process in detail using Monte Carlo and molecu- 
lar dynamics simulations [B[ for a model system, analyzing 
our results in the light of existing theory 0, @, Q. Briefly, 
we discover that a small solid, constrained to remain de- 
fect free by being, at all times, in contact with its own 
liquid (a situation easily realised using optical trapsQ), 
responds to stress by exchanging surface atomic layers 
with the adjacent liquid. Impacting the solid with a mo- 
mentum pulse Q of sufficient strength dislodges an entire 
crystalline layer coherently, which travels into the liquid 
as a distinct though short-lived entity, with a lifetime de- 
termined by the fluid viscosity Q- A curious feature of 
this process is that weaker pulses do not dislodge par- 
tial layers, leading to a novel resonance phenomenon dis- 
tinguished by a pronounced peak in sound absorption Q. 
We believe that our work may be useful for undertanding 
the failure behaviour and sound and heat [HI absorption 
properties of nanostructures 11 , 12] • Since coherent scat- 



tering of momentum pulses occur over a narrow window 
of incident energies, this phenomenon may also be used 
as a detector/analyzer for weak acoustic shocks. Free 
standing, cleaved single atomic layersfHJ have recently 
been shown to posses interesting mechanical and elec- 
trical properties. Coherent spallation^ of nanocrystals, 
as discussed here, may be a practical way to produce 
such atomic layers or for making nanowires or nanosur- 



face coatings[14j in the future. 

We create a dislocation-free nanosolid by trapping a col- 
lection of atoms in their crystalline state within a po- 
tential well of depth —\x over a finite region S placed in 
contact with its fluid (Fig. 1). Trapping of atoms such as 
alkali metals and noble gases may be achieved by optical 
and magneto-optical techniques [8J, using laser powers 
ranging from 1 — 10 2 mW. On the other hand, colloidal 
solids[15| may be manipulated using a number of optical 
techniques lg, 17 1 or surface templates [T^]. 



To see that this trapped nanosolid is dislocation- 
free, we study the equilibrium behaviour of a model 
2-dimensional (2d) nanosolid using a Monte Carlo sim- 
ulation in the constant number (N), area (A) and tem- 
perature (T) ensemble with usual Metropolis moves 0,01 
using the Hamiltonian H — £V Vij + J2 i (p(ri), where 
Vij is the 2-body potential and <f> is the trap potential, 
here approximated as a rectangular potential well of 
depth —fj, (Fig. 1). For numerical convenience we choose 
a system of hard disks (Fig. 1); our main results trivially 
extend to particles interacting with any form of repulsive 
potential, or even when the interactions are augmented 
by a short-range attraction, provided we choose y, deeper 
than the depth of the attractive potential. While all 
our simulations are carried out in 2d, our qualitative 
results should extend to three dimensions (we shall use 
the generic words 'layer' and 'surface' to describe the 
1-dimcnsional line of atoms). The equilibrium density 
profile (Fig. lb) is obtained[20( for different \x at fixed 
average density r\. 

The 'phase diagram', a plot of the density difference 
across the interface Arj/ri = (rj s — ill)/ 1 ! versus fi (Fig. 
2), shows a sharp jump at fi « 8. This transition is 
associated with an entire close-packed layer entering S 
thereby increasing the number of solid layers by one. 
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FIG. 1: (a) A 2d crystal of N 3 atoms confined to a central 
region S of area A a = L x x L s by means of an optical trap 
with a potential cj>(r) — —fi for r £ S which increases sharply 
to zero elsewhere over a width 5. The trapped solid of density 
rjs is in contact with its fluid of density r\i < rj s , such that the 
average density in the entire simulation cell of N atoms occu- 
pying an area A = L x x L y is r\ = nN/AA. The atoms interact 
with a hard disk potential 



is r\ = nl\ 

sia 



Vij — for \vij | > a and 



sets the scale of length and Yij = Yj — the relative position 
vector of the particles. We have chosen 5 — <r/4. We set the 
energy scale by ksT. (b) Equilibrium behaviour for differ- 
ent fi at fixed rj = .699 (density at freezing rjf — .706 [3]) 
obtained by Monte Carlo simulations in the constant NAT 
ensemble with periodic boundary conditions in both direc- 
tions (N = 1200 particles occupy an area A = 22.78 x 59.18 
with the solid occupying the central third of the cell of size 
L s — 19.73). The trap depth [i — 6, supports an equilibrium 
solid of density r/ s = .753 in contact with a fluid of density 
r\i = .672. The closest-packed lines of the solid in S are par- 
allel to the solid -fluid interfaces which lie, at all times, along 
the lines where 4>{y) — * 0. The density profile rj(y) coarse 
grained over strips of width a (averages taken over 10 3 MC 
configurations each separated by 10 3 MCS) , varies from r\i to 
j) a as we move into S. The horizontal lines are predictions 
of a simple free-volume based theory (see Fig. 2) for rj 3 and 
r\i. (c) Superposition of particle configurations from the MC 
run in (b) showing a solid like order (red : high rf) gradually 
vanishing into the fluid (yellow : low rf) across a well defined 
solid-fluid interface. 



The resulting solid is a triangular lattice with a small 
rectangular distortion Sd(r] s , L s ) The qualitative 

features of this phase diagram may be obtained by a 
simple thermodynamic theory (Fig. 2) with harmonic 
distortions of the solid, ignoring contributions from 
spatial variations of the density. We find that the jump 
in the fractional density difference is sensitive to L s and 
vanishes for large L s or 8, the sharpness of the trap. 

Adiabatically cycling the trap depth fi across the jump 
obtains a sharp rectangular hysteresis loop; this indi- 
cates that 'surface' steps (dislocation pairs) nucleated 
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FIG. 2: Plot of the equilibrium fractional density change 
Arj/ri as a function of /i (points (MC data), thick solid 
line (approximate theory)), showing discontinuous jump at 
fi ~ 8. The MC data are obtained by averaging over 10 3 
configurations each separated by 10 3 MCS, while the sys- 
tem is equilibrated for > 10 7 MCS starting from a uni- 
form fluid with density r\ — .699. The approximate theory 
is based on the assumption that a change in /i produces a 
uniform geometric strain Ed from a reference triangular lat- 
tice with the same number of atomic layers. The geomet- 
ric strain Ed is an oscillatory functional of L s , with an am- 
plitude which decays as 1/L S . The Helmholtz free energy 
of the harmonic solid is then given by f s = /a + ^K& £ di 
where /a and Ka are the free energy and Young's modulus 
respectively, of an undistorted triangular lattice which may 
be obtained from simple free- volume theory 0, [22j]. Minimis- 
ing the total free energy density of the fluid + solid regions, 
f = x [f 3 (ris,L s ) — Arjsfi/ir] + (1 — x)fi(rje) with the constraint 
r\ = xr\ s + (1 — x)r/e, where x is the area fraction occupied 
by 5 and fi(rji) is the free energy of the hard disk fluid [23l] 
produces the jump in Arj(n)/ri. Inset (top) shows a cycle- 
averaged hysteresis loop as fi is cycled at the rate of .2 per 
10 6 MCS. (bottom) A plot of the tensile stress jd against 
strain Ed- The arrows show the behaviour of these quantities 
as ji is increased from the points marked A - D. The corre- 
sponding points in the Ari/rj vs. /i plot is also marked for 
comparison. The state of stress in the solid jumps discon- 
tinuously from tensile to compressive from B — » C due to an 
increase in the number of solid layers by one accomplished 
by incorporating particles from the fluid. This transition is 
reversible and the system relaxes from a state of compression 
to tension by ejecting this layer as /i is decreased. 

in the course of adding (or subtracting) a solid layer, 
have a vanishingly short lifetime. Consistent with this 
we find that the jump in Ari(fi)/r/ vanishes when the 
system is minimised at each /i with a constraint that 
the solid contains a single dislocation pair. Interestingly, 
a dislocation pair forced initially into the bulk, rises to 
the solid- fluid interface due to a gain in strain energy 
where they form surface indentations flanked by kink- 
antikink pairs. The confining potential, which prefers 
a flat interface, may remove these indentations either 
by bending lattice layers or annealing the kink-antikink 
pair incorporating particles from the adjacent fluid. The 
second process always costs less energy and happens 
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quickly. The question is how quickly ? 

To study the lifetime of the kink-antikink pairs (surface 
step), we resort to a molecular dynamics (MD) simula- 
tion using a velocity Verlet algorithm [Bj], with the unit 
of time given by r — y/ ma 2 /kgT , where m(= 1) is 
the mass of the hard disks [HI- Using values of m and 
a typical for atomic systems like Ar or Rb, r ~ lps. 
Starting with an equilibrium configuration at /i = 9.6 
(and UbT = 1) corresponding to a 22-layer solid, we 
create a unit surface step of length I by removing a 
few interfacial atoms and 'quench' across the transition 
to fi = 4.8, where the 21-layer solid is stable. A free 
energy audit involving a bulk free energy gain AFLJ, 
going from a 22 to 21 layered solid (Fig. 2), and an 
elastic energy cost oc log(Z) for creating the step, reveals 
that a surface step is stable only if Z > /* ~ 1/L S . 
For small L s , the critical size I* may therefore exceed 
L x , the total length of the interface. Indeed, we 
observe all steps, save a complete removed layer, get 
annealed by particles from the adjacent fluid over a 
time scale of order r. The solid therefore relieves stress 
only by the loss or gain of an entire lattice layer, since 
all other avenues of stress relief entail higher energy costs. 

This mechanism of stress relaxation via the transfer of 
an entire layer of atoms may be exploited for a variety of 
practical applications, provided we can eject this layer of 
atoms deep into the adjoining fluid and enhance its life- 
time. Highly stressed monoatomic layers tend to disinte- 
grate or curl up 13 as they separate off from the parent 
crystal. It may be possible to bypass this eventuality, if 
the time scale of separation is made much smaller than 
the lifetime of the layer. Can acoustic spallation^ be 
used to cleave atomic layers from a metastable, stressed 
nanocrystal? Imagine, therefore, sending in a sharp laser 
(or ultrasonic) pulse, producing a momentum impulse 
(v y (t — 0) = Vq) over a thin region in y spanning the 
length L x of the cell, which results in a weak acoustic 
shock ^ (corresponding to a laser power « 10 2 mW and 
a pulse duration lps for a typical atomic system). 

The initial momentum pulse travels through the solid 
and emerges at the far end (Fig. 3 (a)-(c)) as a 
broadened Gaussian whose width, A, is a measure of 
absorption of the acoustic energy of the pulse due to 
combined dissipation in the liquid, the solid and at 
the interfaces [3, llfj. For large enough pulse strengths 
Vq, this is accompanied by a coherent ejection of the 
(single) outer layer of atoms into the fluid; such co- 
herent nano spallation involves surface stresses of the 
order k B T/a 2 « 1CT 5 N/cm 2 (Fig. 3(d)). In contrast, 
spallation in bulk solids like steel needs acoustic pres- 
sures m excess of 10 5 N/cm 2 usually available only 
during impulsive loading conditions; the ejected layer is 
a "chunk" of the surface. This difference comes about 
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FIG. 3: (a)-(c) Plot of the absolute value of the momentum 
for molecular dynamics times t — .0007 (a), .2828 (b) 
and 2.8284 (c). The green lines show the position of the solid - 
fluid interfaces. The parameters L s = 19.73, n s = -789 and 
fj, = 4.8. The fit to a Gaussian (blue line) is also shown in (c). 
The initial momentum pulse with strength Vo = 6. is given 
within a narrow strip of size ~ er, just to the left of the solid 
region and the Gaussian fitted (and the width A 2 extracted) 
when the maximum of the pulse reaches a fixed distance of 
44.1 from the source. A reflected pulse can also be seen. To 
reduce interference from the reflected pulse through periodic 
boundary conditions, we increase the fluid regions on either 
side, so that for the MD calculations we have a cell of size 
22.78 x 186.98 comprising 3600 particles, (d) A plot of the 
time development of the Fourier component of the local den- 
sity correlation pc(y,t) obtained by averaging, at each time 
slice t, the sum 2j=i n ex p(— i G • (rj — i"i)) over all particle 
positions n within a strip of width ~ a centered about y and 
spanning the system in x. The wavenumber G = (27r/d)n 
where d — .92 is the distance between crystal lines in the di- 
rection n normal to the fluid-solid interface. The solid (cen- 
tral region with pa(y,t) 7^ 0) ejects a layer (shown by an 
arrow) which subsequently dissolves in the fluid. The curves 
from bottom to top correspond to time slices at intervals of 
At = .07 starting from t = 1.06 (bottom). We have shifted 
each curve upward by .03t/At for clarity. Curves such as in 
(a)-(d) are obtained by averaging over 100 — 300 separate runs 
using different realizations of the initial momentum distribu- 
tion. 



because unlike a bulk system, a strained nanocrystal on 
the verge of a transition from a metastable n + 1 to a 
n layered state readily absorbs kinetic energy from the 
pulse. The fact that surface indentations are unstable 
(see Fig. 4(a)) unless of a size comparable to the length 
of the crystal, L x , ensures that a full atomic layer is 
evicted almost always, leading to coherent absorption 
of the pulse energy. The coherence of this absorption 
mechanism is markedly evident in a plot of A 2 against 
Vq which shows a sharp peak (Fig. 4 (b) ). Among the 
two systems studied by us, viz., a metastable (/i = 4.8) 
and a stable (/i = 9.6) 22 layered solid, the former shows 
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FIG. 4: (a) (color on-line) Configuration snapshot from a por- 
tion of our MD cell showing hard disk atoms (green circles) 
at the solid (bottom)- liquid interface (yellow line) as a weak 
momentum pulse (Vo = 2.) emerges into the liquid. The pulse 
initially ejects a few atoms of the interfacial crystalline layer 
(red circles) of a metastable 22 layered solid at fj, = 4.8. The 
resulting large non- uniform elastic strain evidenced by the 
bending of lattice layers, however, causes these atoms to be 
subsequently pulled back into the solid. Only a stronger pulse 
capable of ejecting a complete lattice layer succeeds in reduc- 
ing the number of solid layers by one leading to overall lower 
elastic energy, (b) Plot of the squared width A 2 of the mo- 
mentum pulse after it emerges from the solid as a function 
of Vo for fi = 4.8 ( O) and 9.6 ( □). The solid line is a guide 
to the eye. The absorption of momentum is largest when the 
available kinetic energy of the pulse exactly matches the po- 
tential energy required to eject a layer. The peak in A 2 (Vb) 
so produced is more prominent for the metastable 22 layered 
solid n = 4.8 than for the stable (fi — 9.6) system showing a 
more coherent momentum transfer in the former case. 

a sharper resonance. The eviction of the atomic layer 
is therefore assisted by the strain induced interlayer 
transition and metastability of the 22 layered solid dis- 
cussed above. Spallation is also facillitated if the atomic 
interactions are anisotropic so that attraction within 
layers is stronger than between layers (eg. graphite and 
layered oxides [13[), for our model, purely repulsive, hard 
disk solid, an effective, intralayer attractive potential of 
mean force is induced by the external potential[4j. 

The spallated solid layer emerges from the solid surface 
into the fluid, and travels a distance close to the mean 
free path; whereupon it disintegrates due to viscous dis- 
sipation (Fig. 3 (d)). The lifetime of the layer is around 
2-3 time units (r) which translates to a few ps for typi- 
cal atomic systems. The lifetime increases with decreas- 
ing viscosity of the surrounding fluid. Using the Enskog 
approximation [25| to the hard disk viscosity, we estimate 
that by lowering the fluid density one may increase the 
lifetime by almost three times. The lifetime enhancement 
is even greater if the fluid in contact is a low density 
gas (when the interparticle potential has an attractive 
part[26j). 
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